G3: Genes, Genomes, Genetics

The Candida Genome Database (CGD; www.candidagenome.org) is both a model organism database and a fungal pathogen database. As a model organism database, CGD stores data for Candida albicans, which serves as a model species both for other Candida spp. and for non-Candida fungi that form biofilms and undergo routine morphogenic switching. As a fungal pathogen database, CGD now hosts locus pages for six species of the best-studied pathogenic fungi in the Candida group. Pathogenic Candida species have become increasingly drug resistant and there is thus a pressing need for research into basic Candida biology, epidemiology, phylogeny, and potential new antifungals, as well as a single location where all of the available data are collected, curated, and made easily searchable. CGD curates the gene-based Candida experimental literature in real time, extracting, organizing and standardizing gene annotations. CGD also links clinical data on disease to relevant Literature Topics to improve searchability for clinical researchers. Because CGD curates the literature for multiple species and most research focuses on aspects related to pathogenicity, we focus our curation efforts on assigning Literature Topic tags, collecting detailed mutant phenotype data, and assigning controlled Gene Ontology terms with accompanying evidence codes. Our Summary pages for each locus include the primary name and all aliases for that locus, a description of the gene and/or gene product, detailed ortholog information with links, a synteny view, a JBrowse window with a visual view of the gene on its chromosome, links to Phenotype, Gene Ontology, Interactions, and Expression pages, as well as sequence information, references cited on the summary page itself, and any locus notes. The database also serves as a community hub, where we link to various types of reference material of relevance to Candida researchers, including colleague information, news, and notice of upcoming meetings. We routinely survey the community to learn how the field is evolving and how needs may have changed. Here we describe CGDs new modern web interface and multiple new tools that have been added in the last 6 months, allowing, among other things, users to better understand the available expression data for a locus and seamlessly switch between species for a given locus.

The Tyr::CreERT2 transgenic mouse lines are essential tools for conditional gene manipulation in melanocytes and are widely used in melanoma research. Two independent lines are in common use: the Bosenberg line and the Larue line. Precise knowledge of transgene integration sites is critical for designing complex genetic crosses, yet the integration sites for these lines have only recently been characterised by whole genome sequencing. Here we report that a Tyr::CreERT2 mouse stock routinely used in BrafCA;Ptenflox melanoma models carries a previously undescribed integration on Chromosome 1, distinct from the previously reported Chromosome 2 integration. Using long-read nanopore sequencing, we mapped the integration to an intergenic locus between Alppl2 and Alpi, revealing a 2,430 bp genomic deletion at the insertion site. We developed position-specific junction PCR and qPCR assays to genotype this allele and confirmed that offspring are born at Mendelian ratios. Given that this stock is associated with elevated spontaneous melanoma penetrance, accurate characterisation of this allele has direct implications for the many laboratories employing this widely distributed melanoma model. SIGNIFICANCEThis study identifies a previously undescribed Tyr::CreERT2 transgene integration on Chromosome 1 in JAX strain 013590, a line widely used across the melanoma research community. Accurate characterisation of this allele has direct practical implications for the many laboratories that rely on this model for studies of melanoma initiation and progression.

Genetic interactions contribute to phenotypic variation, yet most studies characterize them in a small number of reference backgrounds rather than across natural populations. Studying genetic interactions in Caenorhabditis elegans is facilitated by RNAi by feeding, which enables gene-specific knockdowns across many genetic backgrounds without the need to construct individual mutant lines. However, individually screening hundreds of isolates across many RNAi conditions is laborious. Here, we developed a pooled platform that uses whole-genome sequencing to quantify the fitness of hundreds of genetically distinct C. elegans strains competing in a single flask. We leveraged this platform to evaluate RNAi responses in a pooled population of RNAi-sensitive C. elegans strains and identified two classes of QTL that influence RNAi responses. The first class influences responses to knockdowns of specific genes and is exemplified by a QTL that is specific to knockdown of mig-6. The second class influences general RNAi-response across multiple knockdowns. One QTL in this class localized to the right arm of chromosome III and was detected by both GWA and xQTL mapping, which suggests that a common large-effect variant underlies this QTL. We identified a T96K missense change in the dsRNA receptor SID-2 as the underlying causal variant. The 96K allele resides in the extracellular domain of SID-2 and confers increased RNAi sensitivity relative to the 96T allele. We hypothesize that the 96K allele enables greater uptake of ingested lumenal dsRNA through higher-affinity binding to the dsRNA phosphate backbone. The platform we developed enables systematic mapping of genetic interactions across natural populations at a previously impractical scale.

The cultivated dahlia has a complex evolutionary history that is likely the driver for its diverse floral forms. Previous studies have proposed cultivated dahlia as an allotetraploid, autotetraploid, and in one study an auto-allooctoploid. Here, we used PacBio HiFi long-reads and Dovetail Omni-C reads to assemble and annotate a genome for a tetraploid Dahlia variabilis cultivar Edna C. We identified Edna C as a segmental tetraploid with signature of auto- and allo-polyploidy events based on subgenome-specific repeat partitioning and single-copy gene tree construction within homoeologous groups. Additionally, we identified a cultivated dahlia-specific centromeric monomer shared between Dahlia variabilis Edna C and the Dahlia pinnata Kelvin Floodlight. Based on finding three major chromosome clusters, we propose that Dahlia variabilis Edna C is a segmental tetraploid with evidence of auto- and allopolyploidy.

During meiosis, accurate chromosome segregation requires significant condensation and compaction. These processes are mediated by condensins, cohesins, and histone tail modifications. We identified that MET-2, a histone methyltransferase that catalyzes the dimethylation of histone H3 lysine 9 (H3K9me2), differentially impacts chromosome size in the male vs. female C. elegans germline. In met-2 null worms, autosomes during spermatogenesis are significantly larger than wild type, while chromosome size during oogenesis is unaffected. X-univalent size in males is also unaffected by loss of MET-2, indicating MET-2 differentially regulates autosomal and X-chromosome compaction in male spermatogenesis. Autosome size is not changed when males harbor a catalytically deficient MET-2 (met-2CD) or have mutations preventing germline histone H3K9 methylation (H3K9R). In addition, met-2 males, in contrast to met-2CD or H3K9R males, have more active RNA pol II in later stages of meiosis. These data suggest MET-2 plays a noncatalytic role in mediating chromosome structure and transcription. In met-2 male germ lines, genes on the X chromosome, which is typically enriched in H3K9me2, are significantly more likely to be upregulated than genes on autosomes, even though X-univalent size is unchanged. These results suggest that MET-2 plays a sex-specific role that is not limited to its enzymatic activity.

Targeted protein turnover of the transcription factor C/EBP is enhanced by the Tribbles (Trib) pseudokinase family of adaptor proteins, which bind to a conserved degron motif in C/EBP substrate and recruit a cognate COP1 E3 ligase to ubiquitinate C/EBP for degradation. We identified a conserved degron in the Drosophila C/EBP ortholog Slbo and fused this motif to a fluorescent protein (GFP or RFP) under the control of a constitutive promoter. Transgenic flies expressing this fluorescent protein (FP)::degron transgene, which we call the Fluorescent Ubiquitination-based Trbl Turnover Indicator (or FUTTI) sensor, showed reduced levels of fluorescence in response to increased levels of Drosophila Tribbles (Trbl) and conversely, increased levels when Trbl activity was reduced. Mutations in the Trbl conserved pseudokinase substrate binding pocket and in a divergent PEST-like motif in the C-terminal tail region disrupted Trbl activity in the FUTTI assay. These data indicate that FUTTI will be a valuable tool to identify sites of Trbl activity in vivo and the critical features of the protein required for function.

Many RNA-binding proteins play pleiotropic roles in the cell, often through disparate domains that can be decoupled. Here we undertook a mutational analysis of yeast RNA-binding protein Ssd1 to better understand domains and processes required for tolerating extra chromosomes, a special class of aneuploidy. Ssd1 is important for cellular fitness when cells carry extra chromosomes, for reasons that are not clear. Ssd1 has multiple distinct domains, many of which are poorly characterized, including multiple RNA binding folds, nuclear import and export signals, an intrinsically disordered domain, an RNA polymerase II-interacting region, and many phosphorylation sites. We measured multiple phenotypes, including the ability to tolerate an extra copy of yeast chromosome XII as a representative aneuploid. Mutating RNA binding domains disrupts fitness in aneuploid, but not euploid, cells and ablates binding of a representative Ssd1 target. However, mutations affecting phase separation or interaction with the kinase Cbk1, which regulates Ssd1 granule formation during the cell cycle, did not affect aneuploidy tolerance. Evaluating correlations across these and other growth phenotypes informed on functional relationships among domains. We discuss implications for Ssd1 function and its role in aneuploidy tolerance in yeast.

The fungal pathogen Colletotrichum sublineola is a major constraint to sorghum production. The sorghum genotype IS18760 is resistant to multiple C. sublineola strains, but the genetic basis of this resistance is unknown. Bulk segregant analysis using whole-genome resequencing of a biparental mapping population identified a major resistance locus, designated ANTHRACNOSE RESISTANCE GENE LOCUS 3 (ARG3), associated with extensive cell death. Fine mapping delimited the ARG3 locus to a 30 kb interval containing three predicted protein-coding genes and one unannotated non-coding RNA gene. Among these, only the non-coding RNA gene showed consistent sequence polymorphisms between resistant and susceptible lines and elevated pathogen-induced expression in the resistant parent, suggesting that it underlies ARG3-mediated resistance. Regions syntenic to ARG3 across grass species contain conserved non-coding sequences enriched for H3K27me3 and H2A.Z chromatin marks, suggesting a conserved regulatory function. The resistance allele is absent from a 111-line sorghum pangenome. The resistant sorghum line harboring ARG3 and the corresponding avirulent pathogen strain originate from the same geographic region, suggesting co-evolution between host resistance and a prevalent pathogen strain. Phylogenetic analysis suggests that loss of resistance in susceptible lines may relate to insertion of a 12.4 kb repetitive sequence. Importantly, silencing of the ARG3 transcript in the resistant parental line abrogated resistance to C. sublineola, confirming that ARG3 non-coding RNA gene underlies the fungal resistance mediated by the ARG3 locus in IS18760. These lines of evidence show that ARG3 is an important target for resistance breeding and provide new insight into non-coding RNA-mediated disease resistance.

Heliococcus summervillei is an emerging invasive mealybug that causes severe dieback in grasses in pastures and turfgrass landscapes. It is widespread in Australia and has recently been detected across the Caribbean, Mexico, and the United States. Accurate identification of mealybugs is challenging due to cryptic morphology, overlapping diagnostic characters, and limited taxonomic expertise and literature, which makes molecular tools essential for regulatory diagnostics and management. We developed the first Cytochrome Oxidase I (COI) barcode for H. summervillei and used it to examine mitochondrial variation across available populations. COI sequences reveal approximately a 10.2% mitochondrial split between the Type A and Type B variants. Phylogenetic, haplotype network, and genetic distance analyses show that all invasive range populations share one haplotype associated with a recent invasion in the United States, Australia, Pakistan, and the Caribbean, whereas the Barbados lineage contains two closely related haplotypes that represent a historically stable mitochondrial variant. Together, these results establish the first COI reference library for H. summervillei, clarify mitochondrial lineage structure, and provide a practical barcode tool that enables rapid identification of invasive populations and supports timely regulatory and pest management responses. Recognizing mitochondrial variants also establishes a framework for resolving lineage-specific biological and management traits and strengthens reconstruction of introduction pathways central to regulatory decision-making and limiting further spread.

Degenerate sex chromosomes (e.g., the Y or W) remain among the most difficult regions of eukaryotic genomes to assemble because they are highly repetitive and structurally complex. While many recent lepidopteran genome assemblies contain W chromosome scaffolds, the accuracy and consistency of these assemblies remain uncertain, due to lack of replication within species. However, the silkworm moth Bombyx mori is an exception, with numerous independent W chromosome assemblies currently available. We compared six independent long-read W chromosome assemblies, which proved to be highly inconsistent in structure, even among nominally identical genotypes. In contrast, autosomes and the Z chromosome were highly concordant among these assemblies, indicating that current assemblies remain unreliable for resolving W chromosome structure. Additionally, we analyzed repetitive DNA content across the genome. First, we combined assembly- and read-based repeat-discovery methods to generate a comprehensive and curated Bombyx repeat library, which we make publicly available. Assessing repeat content and diversity, we find that the W chromosome is comprised almost entirely of repetitive DNA but that the richness and divergence of W repeats are substantially reduced compared to the remainder of the genome. This reduced diversity, initially inferred from assemblies, is confirmed by direct analysis of PacBio HiFi sequencing reads partitioned by chromosome. We also demonstrate that the B. mori p50ma genome assembly (the current NCBI RefSeq assembly) carries a W chromosome and mitochondrial genome introgressed from B. mandarina. This discovery provided an opportunity to investigate patterns of divergence between closely related W haplotypes, revealing substantially more rapid turnover of repeat content on the W than elsewhere in the genome. Together, our results show that current W chromosome assemblies, although structurally flawed, nevertheless capture robust biological patterns of repeat diversity and support the hypothesis that rapid repeat turnover, rather than frequent chromosome replacement, may underlie the apparent lack of W chromosome homology across Lepidoptera.

Gene model for the ortholog of Density regulated protein (DENR) in the Apr. 2013 (BCM-HGSC Dpse_3.0/DpseGB3) Genome Assembly (GenBank Accession: GCA_000001765.2) of Drosophila pseudoobscura. This ortholog was characterized as part of a developing dataset to study the evolution of the Insulin/insulin-like growth factor signaling pathway (IIS) across the genus Drosophila using the Genomics Education Partnership gene annotation protocol for Course-based Undergraduate Research Experiences.

Decades of research has uncovered a wealth of mechanistic information about the development of sensory bristles in Drosophila melanogaster. By studying large-effect, often loss-of-function mutations, many genes have been associated with bristle development, morphology, patterning, and number. Equally, the number of bristles present in certain areas of the fly cuticle is a classic quantitative trait, the genetic basis of which has been studied using a range of tools, from artificial selection to QTL (Quantitative Trait Locus) mapping. Such studies have often implicated well-understood bristle development genes as contributing to natural variation in bristle number. Here we contribute to the study of bristle number genetic variation in flies by executing a GWAS (genomewide association study). We generated whole genome sequences for 897 phenotyped male D. melanogaster individuals derived from a wild-derived, but lab-adapted outbred population, revealing - following quality control and filtering - over 780,000 variants with frequencies greater than 5%. Using these data we estimated the SNP (Single Nucleotide Polymorphism) heritability for ABN (abdominal bristle number) and SBN (sternopleural bristle number) as 0.28 and 0.35, respectively. These values indicate that our set of genotyped variants collectively explain a substantial fraction of the variance in phenotype in the mapping panel. Subsequently, genome scans revealed 1085 (ABN) and 211 (SBN) genomewide significant sites, and - due to extensive LD (Linkage Disequilibrium) in our panel - nearly all these sites are clustered into three locations; We find a GWAS hit for ABN in the middle of chromosome 3L, and hits for SBN at the tip of the X chromosome (where several prior mapping studies have resolved QTL for bristle number), and on 2L. Surveying existing studies that identified genes that control bristle number/development, we highlight several candidates that may segregate for causative, functional variants.

We present a chromosome-level genome assembly of a female Norwegian white-throated dipper (Cinclus cinclus cinclus) generated using Oxford Nanopore Technologies (ONT) long reads and Hi-C scaffolding. The assembly comprises two pseudo-haplotypes, hap1 (1186 Mb) and hap2 (1115 Mb), with 96.7% and 94.4% of sequences assigned to chromosome-scale scaffolds, respectively. Both pseudo-haplotypes contain 40 autosomes, with the Z and W sex chromosomes assigned to hap1. Compared with the PacBio HiFi-based C. c. gularis reference assembly bCinCin1.1.pri, which contains 38 autosomes, sequence represented as a single dot-chromosome (chr 36) is resolved into three distinct dot-chromosomes (chr 36, 39, and 40), a configuration supported by Hi-C contact patterns. BUSCO completeness was high for hap1 (99.2%) and hap2 (95.0%), with 19,003 and 17,746 predicted protein-coding genes, respectively. Compared with the HiFi-based C. c. gularis reference and HiFi-based assemblies generated from the same individual, the ONT-derived assemblies were substantially less fragmented and recovered more sequence from the smallest chromosomes. Synteny was otherwise largely conserved between subspecies. HiFi depletion increased strongly from macrochromosomes to micro- and dot-chromosomes, and HiFi-depleted regions were enriched for repeats and predicted non-B-DNA-associated features, particularly G-quadruplexes and direct repeats, whereas ONT coverage remained comparatively stable. These results show that conventional genome-wide assembly metrics can obscure substantial differences in the recovery of repeat-rich avian dot-chromosomes and highlight the value of chromosome-aware evaluation and ONT sequencing for recovering these regions.

European white oaks (Quercus section Quercus) are ecologically and economically important forest trees characterized by extensive shared genetic variation and a history of interspecific gene flow. Genomic resources remain uneven across species, limiting comparative analyses and pangenome development. Here, we present haplotype-resolved chromosome-scale genome assemblies and genome annotations for four European white oak species: Quercus robur, Q. petraea, Q. pubescens, and Q. frainetto. The assemblies were generated from PacBio HiFi sequencing data and include both phased haplotypes for each species. Genome sizes range from 779 to 817 Mb and all assemblies are organized into 12 chromosome-scale pseudomolecules with high completeness and contiguity. We additionally provide species-specific repeat annotations, structurally and functionally annotated protein-coding gene sets, and complete organellar genomes. The dataset includes the first reference genomes for Q. pubescens and Q. frainetto, together with newly generated assemblies for Q. robur and Q. petraea produced using a consistent sequencing and analysis workflow. These resources provide a standardized framework for comparative genomics, pangenome construction, genome evolution studies, and investigations of adaptation and introgression across European white oaks.

Daphnia magna is a widely used model organism in ecological genomics, ecotoxicology, and evolutionary biology. While a few genome assemblies of this species of freshwater crustacean are available, disagreements in their chromosomal arrangements point toward scaffolding problems. We present the most contiguous chromosome-level genome assembly of D. magna, the first based on PacBio HiFi sequencing technology. The assembly spans 184.38 Mb and has a scaffold N50 of 12.4 Mb. Eight of ten chromosomal scaffolds consist of mainly two contigs, representing one chromosome arm each. Eleven of 20 telomeric regions are fully resolved. Using a genomic map and comparative genomics, we were able to not only resolve the problems of previous assemblies, but also to demonstrate that some of the scaffolding errors followed a systematic pattern, most likely resulting from mis-joins during Hi-C scaffolding. Additionally, we provide expression profiles for annotated genes derived from different D. magna RNA-seq datasets, enabling rapid assessment of predicted gene expression and identification of potential pseudogenes. The new genomic and transcriptomic resource enables comparative genomic studies on genome architecture, structural variation, and gene family evolution in this important model organism.

Long breeding cycles constrain genetic gain in Scots pine, while mature progeny trials can provide reference populations for genomic selection. We combined phenotypic records from two German trials established in 1990 with dense SNP data. We complemented these data with an offspring-level parentage audit, duplicate filtering, and simulations of breeding strategies. Of 1,986 phenotype-matched genotyped trees, 1,668 formed the parentage-green set and 1,651 remained after exclusion of 17 near-duplicate samples. The final population comprised 1,521 supported controlled-cross offspring and 130 mother-known, father-unknown offspring from 87 progeny labels. Across 14 trait-by-age measurements, GBLUP gave the clearest improvements for diameter and volume from age 20 onwards, whereas height was more mixed. At age 35, ABLUP versus GBLUP heritability was 0.155 versus 0.235 for diameter, 0.265 versus 0.263 for height, and 0.172 versus 0.236 for volume. Leave-progeny-out predictive ability increased from 0.205 to 0.238, 0.249 to 0.265, and 0.203 to 0.239, respectively. SNPscan_breeder simulations compared phenotypic selection, progeny testing, cross-generation genomic selection, and genomic selection with phenotypic thinning under five diversity variants. Progeny testing produced the greatest cumulative gain, but its 45-year cycle reduced annual response. Under the base assumptions, genomic selection with phenotypic thinning gave the highest annual gains for height and fungal resistance, whereas pure genomic selection gave the highest annual diameter gain. Genomic strategies accumulated more kinship than conventional strategies, although a {lambda} = 0.10 kinship penalty improved founder retention with little loss of gain; no mitigation option simultaneously maximised gain, prediction accuracy, and diversity. These results support genomic shortlisting within a field-tested programme with explicit reference-population updating and diversity management.

Most genetic diversity stems from spontaneous mutations, that is, errors in DNA repair or replication. But for dozens of organisms across the tree of life, mutations at specific loci are not spontaneous but developmentally programmed: effectively, some organisms edit their own DNA sequences. This is perhaps most common among pathogens and parasites, many of which use editing to diversify genes that produce important antigens. Plant-parasitic potato cyst nematodes are damaging agricultural pests that establish a lifelong feeding site inside the root of their host plant. We previously observed extensive diversity of rare alleles at HYP1, the most highly expressed gene that encodes a protein secreted by potato cyst nematodes during parasitism. Importantly, HYP1 alleles differ from each other by complex, in-frame rearrangements of short repeated sequence motifs within a single exon. Combining several lines of evidence, we previously hypothesized that potato cyst nematodes use developmentally-programmed mutations, or editing, to diversify HYP1 alleles in the soma. In the current work, we now test this hypothesis. We employ highly accurate long-read DNA sequencing of a simplified genetic system to identify potential rare edited alleles, we use a transgenic yeast system to describe large de novo mutations at HYP1, and we interpret our findings in light of key population genetic parameters as well as the genetic diversity surrounding HYP1 and across the genome. Significance StatementIn parasites, genetic diversity is valuable fuel for the coevolutionary arms race with their hosts. However, in parasitic nematode worms, low genetic diversity is often observed. Several other parasites and pathogens (e.g., single-celled eukaryotes) generate genetic diversity in an unusual way: they edit the DNA sequence of an important gene in their own genome, instead of waiting for rare spontaneous mutations to occur. This has never been observed in a plant parasite, but we previously described a major parasitism gene (HYP1) in plant-parasitic nematodes with such diverse DNA sequences that it looks like it could perhaps be edited. Now, we generate better data for a simplified genetic system and show that much of the observed HYP1 diversity was actually cryptic DNA sequencing errors. Although there is real genetic diversity at HYP1, it is best explained not by editing but by fundamental evolutionary forces that operate across the genome.

Italian ryegrass (Lolium multiflorum) is a key temperate forage species underpinning livestock production in Europe. Genomic resources remain limited by its large (2.2 Gb), repetitive, and highly heterozygous genome. Here, we present a high-quality chromosome-scale genome assembly of the Swiss L. multiflorum ecotype Tremona, collected in 2008 in Ticino, Switzerland, and subsequently incorporated into recurrent breeding cycles in the Swiss breeding program. To address systematic assembly artefacts caused by unresolved haplotypes in our initial PacBio HiFi assembly, we developed ParaLies, a post-assembly tool that identifies and removes artefactual duplications based on sequence divergence while preserving true paralogous gene copies. ParaLies reduced the duplicated BUSCO rate from 16.91% to 6.72% without loss of bona fide genomic content. The resulting assembly has a contig N50 of 15.69 Mb and captures 94% of the expected 2.2-Gb genome size. We further analyzed whole-genome resequencing data from Tremona, additional Swiss ecotypes, and publicly available North American germplasm. Tremona was genetically homogeneous, with no evidence of pronounced recent bottlenecks or substantial within-population structure, and was genetically distinct from the other Swiss ecotypes analyzed. Together, the Tremona genome and ParaLies provide valuable resources for L. multiflorum genomics and breeding and demonstrate a scalable approach for reducing haplotype-induced redundancy in highly heterozygous genomes.

Adult stem cells can survive diverse insults, from ionizing radiation to toxic chemicals. In vertebrates, ABC transporter expression and dye-efflux capacity have been associated with stem-cell protection and multidrug resistance. How compound accumulation differs between adult stem cells and their differentiated daughters within an intact tissue remains poorly understood. Here, we show that adult Drosophila intestinal stem/progenitor cells maintain a distinct low-accumulation state for membrane-permeable compounds and that ABC transporters contribute to this phenotype and to the response to cytotoxic stress. Across nine membrane-permeable probes, stem/progenitor cells accumulated eight at lower levels than differentiated enterocytes. We selected tetramethylrhodamine methyl ester (TMRM), a membrane-permeable fluorescent probe that showed a strong difference in accumulation between stem/progenitor cells and enterocytes, as a quantitative readout for an RNAi screen of all 56 Drosophila ABC-family genes. The screen identified seven candidate transporters whose depletion increased the stem/progenitor-to-enterocyte TMRM fluorescence ratio. Further characterization of one of them, ABCB7, showed that ABCB7-depleted guts contained fewer mitotic cells during cytotoxic drug exposure. These results implicate ABC transporters in lower compound accumulation in Drosophila intestinal stem/progenitor cells and support a role for ABC transporters in their response to cytotoxic challenge.

Insects of the subfamily Triatominae, commonly known as kissing bugs, are obligate blood-feeding vectors of Trypanosoma cruzi, the causative agent of Chagas disease. Rhodnius prolixus is among the most epidemiologically important vectors in Latin America, whereas Triatoma rubida frequently invades homes and is a potential vector in the southern United States and northern Mexico. Triatomines likely evolved from predatory reduviid assassin bugs through a transition from feeding on arthropods associated with vertebrate hosts to feeding directly on vertebrate blood. To investigate the genomic basis of this ecological and dietary shift, we generated highly contiguous, near chromosome-level genome assemblies and structural gene annotations for R. prolixus and T. rubida. The new R. prolixus assembly improves scaffold N50 more than 40-fold over the current reference genome, from 1.1 to 43.9 Mb, while reducing assembly gaps by several orders of magnitude. Both assemblies exceed 97% BUSCO completeness. Comparative analyses with representative hemipteran genomes revealed expansions of gene families associated with chemosensation and metabolism, including detoxification, protein degradation, and digestion, together with signatures of positive selection in genes involved in digestive and sensory functions. These assemblies represent the most contiguous and complete genomic resources available for Triatominae and provide a robust foundation for investigating vector biology, host adaptation, and the evolutionary origins of blood feeding within Reduviidae. Interpretive summaryKissing bugs are insects that are known for feeding on blood. They can spread a disease called Chagas disease because they transmit a parasite called Trypanosoma cruzi. To understand how kissing bugs evolved and which genes facilitate blood feeding of vertebrates, a collaboration between scientists at USDA-ARS, University of Georgia, and University of South Bohemia sequenced the genome of two kissing bugs: Rhodnius prolixus and Triatoma rubida. By comparing the genes with those of other insects in the order Hemiptera, scientists discovered that kissing bugs have more genes involved with detecting environmental chemical stimuli and metabolism as well as positive selection for genes involved with digestion and sensory-related proteins. These genome assemblies will help scientists learn more about how these insects evolved, and this research is important for understanding insect feeding biology which can be used to develop methods to control the kissing bugs and the spread of Chagas disease.